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Abstract 

We use the heavy quark expansion to investigate the width difference Ar# s between the B s 
mass eigenstates. The corrections of (D(AQCD/mb) and 0{m s /mb) to the leading order expres- 
sion in the operator product expansion are derived and estimated to yield a sizable reduction 
of the leading result for Ar^ by typically 30%. For completeness we also quantify small 
effects due to penguin operators and CKM suppressed contributions. Based on our results we 
discuss the prediction for {/S.Y /Y) Ba with particular emphasis on theoretical uncertainties. We 
find (Ar/r)s a = 0.16to;Q9, where the large error is dominated by the uncertainty in hadronic 
matrix elements. An accuracy of about 10% in {/ST/T)b s should be within reach, assuming 
continuing progress in lattice calculations. In addition we address phenomenological issues 
and implications of a Ar# s measurement for constraints on AM# s and CKM parameters. We 
further consider in some detail the lifetime ratio t(B s )/t(BA and estimate that, most likely, 
\r(B s )/r(B d ) - 1| < 1%. 
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1 Introduction 



Mixing phenomena in neutral B meson systems provide an important testing ground for 
standard model flavordynamics. The mass difference between the Bd eigenstates, AMs d , gave 
the first evidence for a large top quark mass and provides a valuable constraint on \ V t d\ and 
the CKM unitarity triangle. A direct measurement of AMb s , the corresponding quantity for 
B s mesons, through B s -B s oscillations, would yield further information and help to reduce 
hadronic uncertainties in the extraction of CKM parameters. Complementary insight can be 
gained from the width difference Ar^ between the B s mass eigenstates |L[ [J. This width 
difference is expected to be the largest among bottom hadrons [[J , and it may be large enough 
to be accessible by experiment in the near future. The width difference for Bd mesons, on the 
other hand, is CKM suppressed and experimentally much harder to determine. 

If Ar b s is indeed found to be sizable, the observation of CP violation and the extraction 
of CKM phases from untagged B s data samples can be contemplated Jl], [3], |!|. This possibility 
could be important in two respects. First, tagging any B s data sample costs in statistics and 
in purity. Second, the rapid oscillations dependent on AMs s t all cancel in time evolutions 
of untagged B s data samples, which are governed by the two exponentials exp(— T^t) and 
exp(— Tut) alone. 

The present article continues previous work of one of us [|IJ on the phenomenological 
potential of Ar# s , and focuses on theoretical uncertainties and improvements of the prediction. 
We compute the width difference in the heavy quark expansion and include explicit 1/rrib- 
corrections, which improves over previous estimates of ATb s based on a partonic [§, 0, f|, H |10| 
or exclusive |TI|] approach and allows us to assess the remaining uncertainties more reliably. 
Combined with future measurements of ATb s these predictions can be used to derive indirect 
constraints on |Vt s /V^| @ and AMb s - Non-standard model sources of CP violation in the 
B s system would reduce ATb s compared to its standard model value, as explained in [fT^] , so 
that a lower bound on the standard model prediction is especially interesting. 

Starting from the flavor eigenstates {\B S ), \B S )}, B s — B s mixing is determined by the 2x2 
matrix 

M = M- % -T. (1) 

with hermitian M and T. Due to CPT conservation M\\ = M22 = Mb s , Tii = T22 = r_e s - 
We recall that for the B s system the off-diagonal elements obey the pattern 



1.2 



Mi 



1.2 




(2) 



The mass and lifetime difference between eigenstates are given by ('FT for 'heavy', 'L' for 
'light') 

AM Bs = M H - M L = 2\M 12 \, (3) 

Ar -r r - 2Re ( M i^) U) 



1 



The corrections to (]3|) and (|) are extremely suppressed. They enter only at order |ri 2 /Mi 2 | 2 
and vanish altogether in the limit of exact CP symmetry. Anticipating the actual hierarchy 
of eigenvalues, we have defined both AMg s and Ar# s to be positive quantities. 

Neglecting CP violating corrections, which are very small in the standard model (SM), the 
mass eigenstates are CP eigenstates (up to corrections of at most 10 -3 ), and with the phase 
convention CP\B S ) = -\B S ) one has \B H/L ) = (\B S ) ±\B s ))/y/2. Thenfl, using standard CKM 



phase conventions |13 | 



AT Bs = -2Ti 2 = -2T 21 . (5) 

Note that the lighter state is CP even |I] and decays more rapidly than the heavier state. This 
also follows from the fact that most of the decay products in the b —>■ ccs transition which are 



common to B* and B« are CP even 11 



Both the mass and lifetime difference are determined by the familiar box diagrams that give 
rise to an effective AB = 2 Hamiltonian ('5' denotes 6-quark number). On distance scales 
larger than 1/Myp, but still smaller than l/m&, this effective Hamiltonian contains a local 
AB = 2 interaction as well as a bilocal part constructed from two (local) AB = 1 transitions. 
The mass difference is given by the real part of the box diagram and is dominated by the top 
quark contribution. For this reason, M12 is generated by an interaction that is local already 
on scales x > 1/Mw and theoretically well under control. The short-distance contribution 
has been calculated to next-to- leading order in QCD [fL4| . The long-distance contribution is 



parametrized by the matrix element of a single four-quark operator between B s and B s states. 
Corrections to this result are suppressed by powers of ml/M^ and completely irrelevant for 
all practical purposes. 

The lifetime difference is given by the imaginary part of the box diagram and determined 
by real intermediate states, which correspond to common decay products of B s and B s , so that 
only the bilocal part of the AB = 2 Hamiltonian can contribute. The presence of long-lived 
(on hadronic scales) intermediate states would normally preclude a short-distance treatment 
of the lifetime difference as indeed it does for neutral kaons. But for bottom mesons, the 
b quark mass m;, provides an additional short-distance scale that leads to a large energy 
release (compared to Aq C d) info the intermediate states. Thus, at typical hadronic distances 
x > 1/rrib, the decay is again a local process. The bilocal AB = 2 Hamiltonian can be 
expanded in inverse powers of the heavy quark mass, schematically: 

Imi [ d 4 xT(o AB=1 (x)O AB=1 (0)) = £%(^ B=2 (0) (6) 

The matrix elements of local AB = 2 operators that appear here and in the mass difference 
are not independent of mj,. Their mass dependence could be made explicit with the help of 
Heavy Quark Effective Theory. The difference between the mass and lifetime difference is that 
for the lifetime difference explicit l/m^ corrections arise from the expansion (|6]) even before 
expanding the matrix elements of local operators. The heavy quark expansion applies as well 



t Subsequently, we present the result of our calculation of I^i as a result for Ar^ using (||). If one 
does not want to assume standard model CP violation, (||) must be generalized to (fj), but our result for T21 
is still valid, provided non-standard model CP violation modifies only M12, but not T^. Since Ti2 results 
predominantly from tree decays, this is reasonable to assume. 
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to the diagonal elements Fa = F Bs = (F B + Fl)/2 and has been used to predict the total 
width of bottom hadrons [[|. A contribution to Ti 2 requires that the spectator strange quark 
and the bottom quark come together within a distance 1/rrib in a meson of size 1/Aqcd- This 
volume suppression together with the phase space enhancement, leads to the estimate 



12 



rn 



l 67r *('^£) . (7) 
rrtb 1 



The application of heavy quark expansions to non-leptonic decays assumes local duality. The 
accuracy of this assumption can not be quantified within the framework itself, at least not 
to finite order in the heavy quark expansion. The assumption that the sum over exclusive 
modes is accurately described by the heavy quark expansion might be especially troubling 
for Ar^, since it is saturated by only a few £)(*>**)£)(*>**) intermediate states and the energy 
release is only slightly larger than one GeV. On the other hand, in the small-velocity limit 
Aqcd <^ rrib — 2m c <C m c , and the N c — > 00-limitQ, local duality with only a few intermediate 
states can indeed be verified explicitly [fLlp . 

This article starts from the hypothesis that duality violations should be less than 10% for 
Ar^. Aiming at an accuracy of 10%, the following corrections to the leading order result 
have to be considered: 

(i) 1/mj, corrections from dimension seven operators in (^|). 

(ii) Deviations from the 'vacuum insertion' ('factorization') assumption for matrix elements 
of four-fermion operators. 

(iii) Radiative corrections of order a s /n. 

(iv) Penguin and Cabibbo-suppressed contributions. 

The major part of this paper is devoted to l/m^ corrections. We hope to return to radiative 
corrections in a subsequent publication. These would bring the short-distance part of the 
calculation for AF Bs on the same level that has already been achieved for AM^. The result 
for AF Bs to next-to-leading order in the l/m b expansion is obtained in Sect. |]. We use the 
vacuum insertion approximation for the dimension seven operators, and express the result in 
terms of two non-perturbative parameters that have to be computed with lattice methods. 
Sect. |]is devoted to the phenomenology of Ar^. Numerical results are discussed in Sect. ^ 



together with the theoretical uncertainties in AT Bs /T Bs . In Sect. [372] a generally valid upper 
bound on AF Ba is derived. Sect. |3.3| describes potential strategies to measure the width 
difference in experiment. Some phenomeno logical applications of such a measurement are 
considered in Sect. [T3j. 

An issue related to AF Bs concerns the total decay rate F Bs of B s mesons, averaged over the 
long-lived and short-lived component. For experimental investigations of AF Bs |T| it would be 
helpful to know to what extent the average B s decay rate F Ba differs from F Bd . These decay 
widths are estimated to coincide to a high accuracy ||. We quantify this expectation and 
detail the contributions that could give rise to a difference between F Ba and T Bd in Sect. £|. A 



■^This limit is necessary to justify the factorization assumption for four-fermion operators. 
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Figure 1: Diagram that gives the leading and next-to-leading in 1/mb terms in the heavy 
quark expansion of the forward scattering amplitude. 



summary is presented in Sect. |5|. Penguin and Cabibbo-suppressed contributions turn out to 
shift Ar# s by less then 10% and are discussed in the Appendices, along with a comment on 
the lifetime ratio of B + to mesons. 



2 Ar 



B, 



Basic Formalism 



The optical theorem relates the total decay width of a particle to its forward scattering 
amplitude. The off-diagonal element T 2 i of the decay width matrix is given by 



21 



2M 



■{B S \T\B S ). 



Bs 



The normalization of states is (B S \B S ) = 2EV (conventional relativistic normalization) and 
the transition operator T is defined by 



T = Imi J d A x TH ef f{x)H ef f{0). 



(9) 



Here 7~L e ff is the low energy effective weak Hamiltonian mediating bottom quark decay. The 
component that is relevant for T 2 i reads explicitly 



H 



eff 



7=2 



V* b V cs (C 1 (//)(6 i c i )v_A(cjSi)y- J 4 + C 2 {lj)(biCi)v-A{cjSj)v- 



(10) 



where we are neglecting Cabibbo suppressed channels and the contributions from penguin 
operators, whose coefficients are small numerically. These contributions will be considered 
in the Appendices. We use the notation (q\q2)v-A = <?i7^(l — 75)^2 an d similar notation for 
other combinations of Dirac matrices. The indices i, j refer to color. The Wilson coefficient 
functions C\ 2 read in the leading logarithmic approximation 



C2J 



C + ±C. 



CM 



a s {M ] 



w ) 



-1 6/23 



a s {M ] 



w , 



-1 -12/23 



a s (n) 



(11) 



with scale fi of order m^. 

The leading contribution to the AB = 2 transition operator is shown in Fig. [I], where 
the vertices correspond to the interaction terms in (|i~0p. The operator product expansion is 
constructed using standard methods [§] . Because of the large momentum flowing through the 
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fermion loop, it can be contracted to a point. To leading order in l/m&, the strange momentum 
can be neglected and the b quark momentum identified with the meson momentum. The result 
can be expressed in terms of two dimension six operators 



Qs = {biSi)s-p{bjSj) S -p- 



(12) 
(13) 



The first operator coincides with the single operator that contributes to the mass difference. 
The appearance of a second operator can be traced to the fact that in the calculation of T 2 i 
the external b momentum can not be neglected, because its zero component (in the meson 
rest frame) provides the large momentum scale. 

To include l/m& corrections, the forward scattering amplitude, evaluated between on-shell 
quark states, is expanded in the small strange quark momentum and matched onto operators 
with derivatives or with a factor of m s , the strange quark mass, which we count as Aqcd- 
Operators with additional gluon fields contribute only to corrections of order (Agc^/m^) 2 
and need not be considered. It is more direct (and rather trivial at this order) to use the 



background field method |15| . Since we do not scale out the 'kinematic' part of order m& in 
derivatives acting on b fields, we do not have immediate power counting. Some operators of 
higher dimension in (^) have to be kept, if they contain derivatives on b fields, such as R 2 
below. Using the equations of motion, we are left with operators with at most one derivative 
on b fields and obtain 



21 



G 2 F m\ 



12vr(2M Bs 



-(v; b v cs ) 2 VT^ 



1 - z)K x + -(1 - Az)K 2 ) (Q) + (1 + 2z) (K x - K 2 ) (Q s ) + S x/m 



where z = m 2 /m 2 and 



(14) 



K x = N c Cl + 2C X C 2 K 2 = Cl (15) 

The brackets denote the matrix element of an operator O between a B s and B s state, (O) = 
(B s \0\B s ). The 1/rrif, corrections are summarized in 



$l/r, 



[l + 2z) 
12z 2 



K x {-2{R X ) - 2(R 2 )) + K 2 ((Ro) - 2(R 1 ) - 2(R 2 )) 
K x ((R 2 ) + 2<i? 3 » + K 2 ((R 2 ) + 2(R 3 )) 



l-Az 

The subdominant operators are denoted by Ri and Ri and read (R4 will be needed below) 



(16) 



^0 = Qs + Qs + 7,Q Qs = (biS^s-pfas^s-p 

R\ = —(biSi) S - P (bjSj)s+p 



(17) 
(18) 
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R 2 = -l(6^ p7 "(l - 75)1^(6^(1- lh ) 8j ) (19) 
m h 

R 3 = -1(6^(1 - 75 )D' Si )(6i(l- 7B ) 5j ) (20) 

i? 4 = -1(6,(1 - 75 )^ Si )(6^(l- 75 ) Si ). (21) 

The i?j denote the color-rearranged operators that follow from the expressions for R{ by 
interchanging Sj and s.,-. In deriving (0) we omitted total derivative terms, because four- 
momentum is conserved in the forward scattering amplitude. 

The operators Ri and Ri are not all independent at order l/m&. Relations can be derived 
by using the equations of motion and omitting total derivatives. To reduce Rq, one can start 
from the Fierz identity 



(6 i7M (l - 75)^(6,-7,(1 - 75>i) = (22) 
-(6i7 M (l - 75)5^(^7^(1 - 75 )s f ) + -Qp, (6i7 A (l - 75)sj)(6j7a(1 - 75>i) 
and apply derivatives in an appropriate way. Up to corrections of l/m& (or less), we find 



Rq = 2i?i — i?2 + 2i?4 

-Ro — -Ro 

R 2 = -R2 (23) 
-R3 = R3 + R2/2 

i?4 = i?4 + R\ — R\ — i?2- 

The first of these relations shows explicitly that the matrix element of Rq is 1/ m& suppressed 
compared to Q, which is not directly evident from its definition above. 

At this point, we have expressed the l/m& corrections to Ar# s in terms of five new unknown 
parameters, in addition to the two non-perturbative parameters that appear already at leading 
order, and which also contain implicit 1 / nib corrections. In principle they can all be obtained 
within the framework of lattice gauge theory^. Unfortunately, results accurate to 10% are 
not yet available, especially not for (Qs) (and all the subleading operators). We therefore 
adopt the following strategy: we parametrize the two operators that appear at leading order. 
They can be estimated in vacuum insertion or the large N c limit, but should ultimately be 
computed on the lattice. The operators Ri, Ri, on the other hand, are only of subleading 
importance and we shall content ourselves here with the factorization approximation. 

Following standard conventions we express the matrix elements of Q and Qs in terms of 
the corresponding 'bag' parameters B and B$ 

§ The matrix elements of the subleading operators could be evaluated in the static limit. However, to 
consistently include all 1/mb corrections, (Q) and (Qs) must be computed cither in full QCD or in Heavy Quark 
Effective Theory including 1/mt corrections to the Lagrangian as well as to the effective theory operators. 
The parametrization of 1 /nib corrections to (Q) has been analyzed in ]l6[ . 
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(Qs) 



(Q) 




(25) 



(24) 



where M Bs and fs a are the mass and decay constant of the B s meson and N c is the number of 
colors. The parameters B and B s are defined such that B = B s = 1 corresponds to the factor- 
ization (or 'vacuum insertion') approach, which can provide a first estimate. Factorization of 
four-fermion operators is a controlled approximation only for large N c or for a non-relativistic 
system. In the large N c limit, 5 = 3/4 and Bs = 6/5. In the sense of these limiting cases, 
factorization for realistic B s mesons can be expected to yield the correct order of magnitude 
and, in particular, the right sign of these matrix elements. Existing nonperturbative calcu- 
lations like lattice simulations for (Q) , and for its counterpart in the K — K system, are in 
agreement with this expectation. Beyond these limits factorization does not reproduce the 
correct renormalization scale and scheme dependence, necessary to cancel the corresponding, 
unphysical dependences in the Wilson coefficients. This raises the additional question, to 
which we return below, at what scale factorization should be employed to estimate the matrix 
elements. Without further information a certain variation of the parameters B, B s should be 
allowed in performing a numerical analysis. 

Next we consider the subleading operators Ri, Ri, where we apply factorization. Using 
relations such as (a, (5 refers to spinor indices, i,j to color as before) 



valid to first order in l/m&, all matrix elements can be expressed in terms of fs s , Mb s and 
quark masses. We find 




1 



(26) 






(27) 






7 



Combining the above results, one can obtain ATb s from fll^p. The sensitivity to V c b may 
be eliminated by normalizing to the total decay rate Tb 3 expressed in terms of the semileptonic 
width and branching ratio 



T(B S -> Xev) _ G\m\ 2 g(z)n QCD 



V A \' ?1~' '"^ r , (28) 



B(B S -> Xez/) 192vr 3 1 1 5(5, -> Xev) 

g(z) = l-8z + 8z 3 - -12z 2 \nz, (29) 

where B(B S — > Xev) is to be taken from experiment^] and z = m 2 c /m 2 as before. t\qcd denotes 
the one- loop QCD corrections (m& refers to the 6-quark pole mass). Their analytic expression 
can be found in fl7fl . At m h = 4.8 GeV, m c = 1.4 GeV, /x = m&, and with a s {m b ) = 0.216 
one has t/qcd = 0.88. Since radiative corrections to ATb s are not yet known, the inclusion 
of radiative corrections to the semileptonic width seems somewhat arbitrary. On the other 
hand, with V c b = 0.04 and r# = 1.54 ps, one obtains m;, »s 4.8 GeV from (|28|) , compared to 
nib ~ 4.5 GeV without QCD corrections. We prefer the first value as our central choice for m& 
in the numerical analysis, but repeat that, in the absence of radiative corrections to ATb s , 
fjQCD can as well be considered as a normalization uncertainty that replaces the normalization 
uncertainty due to the errors in V c b and r# s . Finally one arrives at the following expression: 



— ^ = 16tc 2 B (B s - Xev ^— • 

Tb s 'g{z)r lQCD ml 

(2(1 - z)K x + (1 - 4z)tf 2 ) (l + i-) 5 (30) 

Ml /„ 1 



+ (1 + 2.) (* a - ^2 - W JB 8 + S 1/m + ^ 

5x/ m is related to <5i/ m , defined in flTED , through 

Sl/m = fB 8 M BM/rn, (31) 

and from now on we imply that Q2"7| ) is used. We have indicated by S rem the contributions from 
CKM-suppressed intermediate states (uc, uc, uu) and from penguin operators in the AB = 1 
effective Hamiltonian, which are estimated in the Appendices A and B to be below ±3% and 
about —5%, respectively, relative to the leading order contribution. We shall neglect 5 rem in 
the analysis to follow. 

Since f Bs ~ A 3 ^ D /ml /2 , AT Bs /T Bs ~ 167r 2 (A QC - D /m fe ) 3 as in the estimate (0). Eq. (|30|) is 
valid to leading ((9(l/m 3 )) and next-to-leading order (0(l/mf)) in the heavy quark expansion. 
The most important neglected terms are radiative corrections of order 0(a s /m b y ). Implicit 
here is the assumption that the quantity (AT/T)b s can indeed be represented to reasonable 
accuracy by the series in powers of AQc^/m;, that is generated by the heavy quark expansion. 
As mentioned earlier, this assumption is equivalent to the assumption of local quark hadron 
duality. 

^ Since we show in Sect. ^ that the lifetime difference between B s and Bd is tiny, no attention has to be 
paid to the flavor content of the B meson. 
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The leading term in fl3"0|), represented by the contributions proportional to B and Bs, 
agrees with the results that have been given previously in the literature^ g |7|, §, g, Note 
that we have consistently kept the distinction between quark masses, arising from the short- 
distance loops or the equations of motion, and the meson mass Mb s from hadronic matrix 
elements, since we are aiming at effects beyond leading order in the heavy quark expansion. 

In (^), Ki,K 2 and B,B S should be evaluated at a scale of order m b . If we wanted to 
use vacuum insertion to estimate the bag factors, it is physically clear, especially in the heavy 
quark limit nib — > oo, that vacuum insertion should be applied not at the scale m b , but at 
a typical hadronic scale fi^ ~ 1 GeV. This still leaves us with an ambiguity as to the choice 
of \ih and in addition with the question, how B(u>h) = Bs(fJ^h) = 1 are related to B(nib) 
and Bs{mb). This latter question can be answered in the limit /i^ <C m;, and corresponds to 
the inclusion of 'hybrid logarithms' 0J, as done in [flO | . The evolution from nib to fih is 



performed in the leading logarithmic approximation in the static theory and leads tof^] 



B(m b ) = 1 (32) 

8/25N 



B s (m b ) 




a s [m b ) 



The first equation in fl32|) reflects the well-known result that the matrix element of the operator 
Q has the same leading logarithmic corrections in the static theory (HQET) as the square of 
the decay constant, f Bg . Taking jj, h = 0.5, 1,2 GeV results in B s (m b ) = 0.80,0.88,0.94. (The 
scale Hh = 0.5 GeV might already be too low for a perturbative evolution.) 

The 6-quark mass nib k, 4.8 GeV is probably not large enough to make this estimate real- 
istic, even if factorization held at the scale fih- The logarithm Innib/f^h is not very large, so 
that other contributions like non-logarithmic 0(a s ) terms which are omitted in (|3~2"D, can be 
expected to be numerically of the same order as the hybrid logarithms that are retained, espe- 
cially since summing hybrid logarithms amounts to a moderate 10% effect (with = 1 GeV). 
The one-loop matching of Q on its counterpart (s) in Heavy Quark Effective Theory indeed 
exhibits sizeable cancellations between logarithms and constants, at least in the particular 



matching scheme considered in |20]. Furthermore, the QCD renormalization between nib and 
Hh in (^) is only valid at leading order in HQET and neglects 1 / m b corrections in the matrix 
elements, which is not consistent with our keeping of explicit 1 /nib corrections. On the other 
hand the B factors are in principle calculable in full QCD. In this case they will automatically 
include 1 /nib corrections as well as the hybrid logarithms, among further important contri- 
butions. For these reasons we prefer to keep the expression for (AT/T)b s in the form given 
in ( |30"D and do not include hybrid renormalization explicitly, with the understanding that the 
bag factors will eventually be available from lattice QCD. In our numerical analysis, we take 
the conservative, but perhaps too agnostic attitude that Bs{nib) could take any value between 
0.7 and 1.3, keeping in mind (^) as a particular model estimate of B and Bs- The upper end 
of this range is motivated by the N c — > oo limit, in which Bs = 6/5. 

"Often factorization is assumed for the leading order term, so that B and Bs have to be set to unity to 
recover the result. 

**We have checked the calculation of hybrid logarithms and agree with the findings of p0[. 
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m b /GeV 




a 


b 


c 


(Ar/r) Bs 


4.8 


m b 


0.009 


0.211 


-0.065 


0.155 


4.6 


m b 


0.015 


0.239 


-0.096 


0.158 


5.0 


m b 


0.004 


0.187 


-0.039 


0.151 


4.8 


2m b 


0.017 


0.181 


-0.058 


0.140 


4.8 


m b /2 


0.006 


0.251 


-0.076 


0.181 



Table 1: Dependence of a, b and c on the 6-quark mass and renormalization scale for fixed 
values of all other short- distance parameters. The last column gives (AT /T) Be for B = B$ = 1 
(at the given scale fi), fs s = 210 MeV. 



3 Ar# s — Phenomenology 

3.1 Numerical Analysis of (AT/T) Bs 

We first turn to a numerical analysis and discussion of (AT/T)b b based on (|3~0f). It is useful to 
separate the dependence on the long-distance parameters fg s , B and Bs and write (AT /Y)b s 
as 

(*£\ =\nB + bB„ + r](-L 



V r 



B. 



aB + bB s + c — ^ — , (33) 
1 210 MeV/ ' V ; 



where c incorporates the explicit l/m b corrections. To estimate the sensitivity of (AT/T)b s 
on the short- distance input parameters, we keep the following parameters fixed: m b — m c = 
3.4 GeV, m s = 200 MeV, = 200 MeV. In addition M Bs = 5.37 GeV and the semileptonic 
branching ratio is B(B S — ► Xev) = 10.4%. Then a, b and c depend only on m b and the 
renormalization scale fi. For some values of m b and fi, the coefficients a, b, c are listed in 
Table [l]. For a central choice of parameters, which we take as m b = 4.8 GeV, /i = m b , 
B = B s = 1 and f Bs =210 MeV, we obtain 



(£). 



= 0.220 - 0.065 = 0.155, (34) 

where the leading term and the l/m b correction are separately displayed. As seen from 
the Table, the dependence on m b is weak, but (AT/T) Ba increases by almost 20% when 
the renormalization scale is lowered to m b /2, at fixed B and Bs- These dependences are 
not specific to the values B = Bs = 1. The weak m b dependence is a somewhat accidental 
consequence of using the semileptonic branching ratio to eliminate V cb . If instead we normalize 
to T^ 1 = 1.54 ps and take V cb = 0.04, (Ar/r) Ss would vary from 0.143 to 0.166 under the 
same variation of m b as in the Table. Let us also add the following more general observations: 
(i) The theoretical expression for ATb s in fl30D predicts the sign of this quantity, which a 
priori could have either value. Ar^ 3 is positive and implies a larger decay rate for the CP even 



(lighter) state [10, | ![ (see the conventions in the Introduction). The typical magnitude of 



(AT/T)b s to leading order in the heavy quark expansion is about 0.2, larger than other width 
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differences among bottom hadrons with the possible exception of the case of A& (depending 
on whether theory or present experiments turn out to be right on A&). 

(ii) The explicit l/m& corrections are numerically important and vary strongly with m&. 
For our central parameter choice they reduce the leading order prediction by about 30%. 
Essentially all the various 1 / m b correction terms add with the same sign and make the result 
somewhat larger than the natural size of the corrections, h.QCD/ m b ~ (Mb s — vn b — m^)jm b pa 
8% and m s /m b pa 4%. 

(iii) The contribution from the scalar operator Q$ by far dominates over the contribution 
from Q, because there is a strong cancellation between terms of different sign in the Wilson 
coefficient of the latter operator. This has important implications for (AM/ Ar) b s , which we 
discuss below, because hadronic uncertainties cancel only partially in the ratio B/Bs- 

(iv) If Bs = 1.3, a (AT/T)b s of as much as 0.25 is not excluded, although this appears 
unlikely. On the other hand, if Bs < 1, as suggested by the estimate from hybrid logarithms, 
and if /# 3 turns out to be merely 180 MeV, ( y AT/T) Bs could be as small as 0.07, making its 
experimental detection more difficult. 

This discussion shows that to resolve the theoretical uncertainties, a reliable calculation of 
Bs is mandatory. Further improvement then requires a full next-to-leading order calculation 
of short-distance corrections. 



3.2 Upper Limit on AT Bs 

Since the b — ► ccs transition is the dominant contributor to (Ar)# s , one obtains the upper 
bound 0, [2TJ 

/ I API \ 

ccs) Bs ■ (35) 

It can be readily understood by considering the limit in which only b — *> ccs transitions were 
generated by the effective Hamiltonian. Eq. ([35]) then follows from the requirement that the 
decay rates be non-negative, r± = F(b — > ccs) ± Ar/2 > 0. B{b — > ccs) Bs denotes the fraction 
of 5 s -meson decays governed by the b — > ccs transitions in the absence of mixing. CLEO J2"2] 




recently confirmed our prediction |23| of a significant 'wrong' charm yield in B decays, thereby 



completing the first direct measurement of 

B(b -»■ ccs 1 ) pa B(b -> c) = 0.227 ± 0.035, (36) 

where B{b — > c) is the average number of c produced per b decay. The Cabibbo allowed 
transition is 

B(b -> ccs) = \V CS \ 2 B(b -> ccs') = 0.22 ± 0.03, (37) 
Assuming B{b — > ccs) Bs ~ B(b — > ccs) then yields the upper limit 



|Ar| 



| < 0.44 ±0.06. (3? 

B s 

-3 



Within the heavy quark expansion, (|Ar|)/r)s s is suppressed by m b relative to spectator 
branching ratios, such as B(b — > ccs). From this point of view a bound like (53) might appear 
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trivial. However, the virtue of relation ( |3"5| ) is its very general validity. It would hold even if 
a heavy quark expansion were not applicable to the underlying process. 



3.3 Measuring AT Bs 

We hope to have convinced the reader about the importance of an accurate measurement of 
Ar. One method is to substitute Ts d for the average B s width Tb 3 and to extract ATb s from 
the time-dependences of untagged flavor specific B s data samples Q . Time-dependent studies 

of angular distributions of untagged B s ^ J/'4 ,( j> decays allow the extraction of I\, and also of 
Th if the CP-odd component is non- negligible || |24| . These and other methods using decay 
length distributions of fully reconstructed B s mesons are at present statistics limited [I], [S], |2l| . 
As an illustration one may consider the measurement of 



T 



(B s - JIM) 



L34 +o.23 



-0.19 



±0.05ps 



(39) 



(-) 



recently obtained by the CDF collaboration from a single lifetime fit of their B s — > J/ip4> data 
sample f25fl. Next we can write 



l/r L < r{B s 

which holds only as an inequality, because B s — 



- J (40) 
J I i\x$> is not necessarily a pure CP-even final 



state. The world average Bd lifetime [26 



r Bd = 1.54 ± 0.04 ps 



(41) 



together with the result of section |], informs us about the inverse of the average B s width 
1/Tb s = TB d . We then use 



and obtain 



Ar 



2(^-1 



> 0.3 ±0.4, 



(42) 



(43) 



which is still inconclusive, but can serve to indicate the present status. 

Just establishing a non-vanishing difference in decay length distributions for partially re- 
constructed B s mesons in comparison to the other B mesons would constitute progress. The 
ideal inclusive 6-hadron data sample should have large statistics and be highly enriched in B s 
decay products originating predominanty from a single mass eigenstate B^ (or Bu). The last 
requirement maximizes differentiation between B s and other I?-mesons. The <jxj)X final state 
serves as an example \2%. The probable decay chain is B s — ► D+DjX, which is dominantly 
CP even ||11|| . Both -D s 's then decay into 0's. While D s is seen significantly in 0's, the D + is 
seen in </>'s by about a factor of 10 less and the D° even less than that p8|. The background 



due to 5-meson decays is thus controllable and further suppressed because B's prefer to be 
seen as D° over D + by a ratio of 2.7 [p£]. If sufficient statistics is available, the DfipX sample 
would be even better. 
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The inclusive B s — > <p£ + X sample with a high Pr,rei lepton, is flavor specific. Its time 
dependence is governed by the sum of two exponentials, exp (— r^t) + exp (— r#t). Theory- 
predicts (Tl + Th)/2 = t/T~B d , but the observation of the two exponents requires precise decay 
length and boost information, whose accuracy increases the more fully the B s is reconstructed. 

The less reconstructed the B s data sample, the more important it is to have a mono- 
energetic source of B s mesons. Thus the more inclusive techniques tend to be more useful for 
e + e~ — > Z° experiments than at hadron accelerators. Of course, fully reconstructed B s data 
samples allow clean measurements of Ar# s . 



3.4 B s - B s Mixing and CKM Elements 

The traditional methods for observation of CP violation and the extraction of CKM phases 



require to resolve the rapid AM Ba t oscillations of tagged B s data samples P0| . Current 
vertexing technology allows to resolve such oscillations for AMg s ^ 10ps _1 . Thus the recent 
lower limit from the ALEPH collaboration pT 



AM B , >6.6ps _1 (95%C.L.) (44) 

is significant. It may indicate the need to develop new methods capable of higher resolving 
power. Reliable predictions of AM Bs are therefore important in order to plan future B s 
experiments, in particular if only lower limits will be available with current vertex techniques. 
The most straightforward method makes use of [j3^| 



AM, 



B s 



G 2 F M*y 

6tt 2 



TIB So{x t )MB 3 B B J 2 B\V t 



(45) 



where Xt = mf/M^r. The current relative uncertainty is about 50% and is dominated by the 
uncertainty in B Ba (±30%), f Bs (±40%), \V ts \ 2 (±15%) and S (x t ) (±8%). The fractional 
uncertainty on AM Ba can be expected to decrease to ~ 15% by the year 2002, anticipating 
improvements in the accuracy of the relevant parameters B Ba (±10%), f Bg (±5%), \V ts \ 2 
(±5%) and So(x t ) (±3%). 

A variant of this method uses the experimental value for AM^ and the ratio 



(AM) 



(AM) Bd M Bd B Bd f 2 Bd 



V t 



ts 



Vt 



td 



(46) 



to predict AMb s - This approach will be useful only if the CKM ratio |V^/14i| is accurately 
known. 

If the first observation of B s — B s mixing is a nonvanishing Ar# s rather than AMb b , then 
a complementary method to predict AM Bs opens up, based on the quantity (see (|30"D ) 



Ar \ 

AMJbs 



7r mr 



2 M$y 



V cb V cs 



v ts v tb 



VBSo(x t ) 



(2(1 - z)K t + (1 - Az)K 2 
+ (1 + 2z) (K 2 - K x 



1 



(m b + m s 



(47) 



1 \ B s , , 
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This result is valid to next-to- leading order in the 1 / m& expansion and to leading logarithmic 
accuracy in QCD. We have again used factorization for the subleading l/m& corrections. Note 
that with the bag parameter B as defined in (p4|), the appropriate QCD correction factor t)b 
is identical to C + (/i) from ( pTT| ) in the leading logarithmic approximation. 

In the ratio AT /AM the decay constant cancels and the CKM uncertainty is almost 
completely removed since 

l '-''*'" s " 1±0.03. (4* 



V ts V tb 

At present the accuracy of AT / AM is still rather poor, AT / AM = (5.6 ± 2.6) • 10~ 3 . The 
breakdown of errors is as follows: ±2.3 from varying Bs/B between 0.7 and 1.3, ~t ' 7 from 
varying /i between m^jl and 2rrib, ±0.4 from m& = 4.8 ± 0.2 GeV and ±0.4 from m t = 
176 ±9 GeV. The dominant uncertainty is due to Bs/B, which has never been studied before. 
It is conceivable that a lattice study could actually calculate Bs/B more accurately than the 
bag parameters themselves, because some systematic uncertainties may be expected to cancel 
in the ratio. The quantity AT / AM might thus be calculable rather precisely in the future and 
AM Ba could then be estimated from the observed Ar Bs . In conjunction with AM Bd this would 
provide an alternative way of determining the CKM ratio |V^ S /V^|, especially if the latter is 
around its largest currently allowed value ||. The width difference, and hence its observability 
increases the larger \V ts \ ~ \V c i,\ becomes. In contrast, the ratio T(B — > K* r y)/T(B — > {g, 
is best suited for extracting small |Ks/Kd| ratios, provided the long distance effects can be 
sufficiently well understood P5| . 

These approaches could complement other methods to determine |V^/Vt s |. Such additional 
possibilities would be to relate \V ts \ to the accurate |V^| measurements and to obtain \V t d\ 
from AMs d i CKM unitarity constraints |54| , and in particular B(K + — > tt + vv) which 



has the unique advantage of being exceptionally clean from a theoretical point of view. 



4 The B s - B d Width Difference 

The ratio of the B s and Bd meson decay widths Ts s /TB d is expected to be very close to 
unity p, [3(J. Deviations arise predominantly from SU(3) breaking effects in already small 
corrections to the leading spectator decay of the bottom quark. In the following we will 
discuss the mechanisms that differentiate between r# s and and estimate their numerical 
importance. The decay rate of B^, B s mesons has the general form (q = d, s) 

T Bq = T + Aril + AT%l g + ATg> A . (49) 

Here To denotes the leading, universal free 6-quark decay rate, ATkm is the time dilatation 
correction, Ar mag the contribution from the chromomagnetic interaction of the heavy quark 
spin, and Ar^^ describes the weak annihilation of b with q. Ar^ n and Ar ma5 are of the 
order 0(l/ml) relative to T and AIV^ enters at order 0(l/m\). Higher orders have been 
neglected in (f49|). There is no linear correction in 1/rrib @. Through order 0(1/ ml) one may 
thus write 

r Ba 1 | Arjj -Arjj , Arw.-Arw, | Arj&-Ar{& (5Q) 
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We will now discuss the three different corrections which contribute to T^/r^ — 1 in turn. 
The first two can be related to meson mass differences. For this purpose we define 



M H = ~(M H + 3M H *] 



(51) 



where M# and Mff* are the masses of a pseudoscalar heavy-light meson H ( 1 5'o) and of its 
vector meson partner H* ( 3 S\). In the weighted average Mjj the spin splitting contribution 
cancels in the HQET mass formula which then takes the form (Q = b, c) 



M Hq =m Q +A q + 



2m c 



+ o 



^■QCD 

m 2 r, 



(52) 



Here (p 2 )^ is the average momentum squared of the heavy quark inside the meson and A 9 may 
be viewed as the constituent mass of the light degrees of freedom. Both quantities depend on 
the light quark flavor q but are independent of the heavy quark mass. Combining (|5^) for the 
cases of D S1 D + , B s and B d and recalling that AT^ n /T = — (p i ) q /(2m b ) one finds 



ZAi kin LAL kin 



m c /m b 
m b - m c 



M 



M 



(Mi 



M 



B d/ 



-(3±6)-10" 4 . (53) 



All required meson masses can be obtained from [13[], except for Mb*- In this case we use the 
heavy quark symmetry relation 



M B * - M Ba 



Mb* - M 



■(M B * - M B J = (46 ± 1) MeV 



(54) 



M,y - M D+ 

to find Mb* = (5421 ± 6) MeV. This expectation is in accordance with direct measurements 
of the B* ^ 5,7 transition, which yield M B , - M Bs = (47.0 ± 2.6) MeV |37j. We see that the 
correction in fl53"D is exceedingly small. This number, however, should probably not be taken 
at face value. Given the smallness of the effect it is conceivable that terms neglected in ( |52"D 
could have an impact on the precise estimate of (53). The typical size of such a correction 
would be (here we use Aqcd = 0.3 GeV) 



ZAi kin ^ kin 



m c /m b 
m b - m c 



A 3 
ml 



12 • 10" 



(55) 



At any rate, while 



might not be a completely accurate estimate of this correction, it 



seems safe to conclude that the effect on r^/T^ due to AT^ — AT^„ is well below 1% and 
thus negligible for all practical purposes. 

Next, the chromomagnetic correction AT^ ag can be related to the spin splitting in S-wave 
B mesons and is proportional to Mb* — Mb ■ Hence we may write 



mag 



mag 



ArW, M bi - M Bs - (Mb* - M Bd , 



M, 



B* 



Mb. 



-(3 ±8) ■ 10" 



(56) 



The quantity AT$ ag /T is known and can be calculated to be —0.012. Using Mb* — Ms d = 
(46.0 ±0.6) MeV @ and M B * -M Bs = (47.0 ±2.6) MeV one finds the numerical estimate 
quoted in (p6|). Clearly this effect on the B s — B d lifetime difference is negligible as well. 
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Finally, we turn to the corrections due to weak annihilation. These contributions arise 
from the annihilation reactions bs — > cc and bd — > cm in the case of a _B S and a -B^ meson, 
respectively. Neglecting Cabibbo suppressed modes and penguin contributions they are readily 
calculated to beQ 



Ar 



m 2 



-(l-^)(K lJ B; s) + — K 2 S 2 



Vqcd 

— Kn B^ 

N, 



(57) 



+ (1 + 2z) K\B 



1 

iv; 



Ar 



16n 2 B(B -> Xei/ 



/I, 



1 T/ 2 



(5? 



-(1 + 



11 N r 2 



1 

iv; 



Here we have again used ( p8f) to eliminate the dependence. The leading log QCD co- 
efficients K 12 are defined in ([T^). The bag factors i?^ parametrize the following matrix 
elements 



{B q \(biqi) V - A (qjbj)v-A\B q ) 
{B q \^iqj)v-A%h)v-A\B q ) 
{BqKbiq^s-piqjb^s+plBq} 
(B q \(biqj)s-p(qjbi)s+p\B q ) 

where we have assumed M Bq ~ ra^. 

Using the strict factorization estimate B^ 
f Bd ^ f Bg and expanding in z ~ 0.1) 



)B a m bB\ 



(?) 



N, 



flrnlB 



(?) 



(?) 



(59) 



1 would yield the following result (taking 



af (s) - Ar (d) 



24n 2 B(B -> Xez/) 



1 - 2z 



fact. 



™>t 9{z) Vqcd 



(60) 



Note that, in 'vacuum insertion', this expression coincides with AT^/r while Ar[y A /T is 
twice as large. For our central parameter set, eq. (|60"D amounts to 2- 10 -4 . The extreme small- 
ness of this number is the result of two effects. The first is helicity suppression, manifesting 
itself in the factor of z = m 2 /m 2 in (|60[) . Secondly, a further suppression comes from a - some- 
what accidental - cancellation between QCD coefficients in K\ + K 2 /3 ~ —0.39 + 0.42 = 0.03. 
It is important to realize that both features are a consequence of the factorization assumption. 
Even with small deviations from factorization the factor z{K\ + K%jN^) would be substituted 

Our results are in agreement with the expressions recently obtained in |36|. 
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by a number almost one hundred times larger. To get an idea of the typical order of magnitude, 
we approximate (|57j) and ( ^8[) to 



Ar 



(<?) 

WA 



1Qti 2 B(B -> Xeu 



™6 9{z)flQCD 



(61) 



1 



+ ^(^ " #2 



(<?)> 



+ 0(z) 



where we have used y/1 — 4z ~ (1 — z) 2 « 1 — 2z and neglected small helicity suppressed 
contributions proportional to 2 in the square brackets. Taking —K\ ~ K 2 /3 ~ 0.4 and 



l r( 9 ) p 



(a) 1 



< 0.6, the modulus of the term in square brackets is 0.5 or less, 



which yields AT^/ A /V < 0.023. Assuming 40% of SU(3) breaking then gives 



AT^ — AF 
ZAi WA lAL 



(d) 
WA 



< 0.9%. 



(62) 



Although with extreme variations, allowing also \K\\ and |i^2/3| to differ (for example by 
choosing a renormalization scale \x different from m&), this difference could be up to 2.5%, it 
is more likely that the correction fl6l|) will actually be much smaller due to various possible 
cancellations in ( |6T| ) and because 40% is probably an overestimate of the magnitude of SU (3) 
breaking. Furthermore, from previous experience with lattice calculations of bag parameters 
in the B meson system it seems likely that the B^ will not differ too dramatically from one, 
so that ([62]) , although admittedly somewhat crude, is probably on the safe side. 

Summarizing the discussion of the various contributions to (BO) we conclude that, most 
likely, the ratio of rates of B s and Ba mesons should differ from unity by no more than one 
percent 



1 



< 1%. 



(63) 



5 Summary 

In this paper we have analyzed the theoretical prediction for Ar# s within the framework of the 
heavy quark expansion. We have calculated the explicit next-to-leading 0(l/m&) corrections 
in the operator product expansion for the transition matrix element. In addition to the two 
leading dimension-six operators, five new operators of dimension seven appear at this level. 
The matrix elements of the latter operators were evaluated using factorization, which should 
give a fair estimate of these subleading corrections. Their effect on Ar^ 3 , formally of order 
0(AQ C D/ m b) and 0(m s /m b ), turned out to be sizable numerically, causing a 30% reduction 
of the leading order prediction. 

We performed a numerical investigation of Ar b s with emphasis on theoretical errors, which 
are presently dominated by the uncertainties in hadronic matrix elements. These errors are 
still rather large and lead to a prediction of (AT/T) Bs = 0.16^o;og. However, a systematic 
improvement of this result is possible, in particular by progress in lattice QCD. In the future 
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it would be desireable to measure on the lattice the S-P four-fermion operator along with 
the V-A operator that has received most attention in the past due to its connection with 
the mass difference. Eventually an accuracy of 10% for AT Bs should be feasible when the 
next-to-leading analysis of short-distance corrections is also completed. 

The effects of penguin operators and contributions from CKM suppressed modes have also 
been considered. They were shown to give only a few percent relative correction in (AT/T) Bs 
and are thus negligible in view of the other uncertainties. 

We further studied the B s — B d lifetime difference and quantified the expectation r Bs « r Bd , 
estimating \t(B s )/t(B c i) — 1| < 1%. This result is useful input for experimental analyses of 

To put our theoretical analysis into perspective, we have included a short discussion of the 
current experimental situation concerning AT Bg . Using information on r(B s — > J/i(j<p) and 
t(B s ) = t(B^), we have attempted a preliminary extraction of AT Bs , obtaining (AT/T) Bs > 
0.3 ± 0.4. This is still inconclusive but can be improved by better statistics in the future. We 
have also proposed an alternative route towards a measurement of AT Bg that makes use of 
the (jxpX and/or Df<pX final states in B s decay, which are expected to be dominantly CP 
even. The present experimental information may be complemented by the bound (AT/T) Bs < 
2B(b — > ccs)b s ~ 0.44 ±0.06. This bound is not very strong, but it has the advantage of being 
valid independently of the heavy quark expansion and it is interesting for principal reasons. 

In addition we have briefly reviewed some phenomenological applications that could be 
opened up by further progress on the experimental as well as the theoretical side. These 
possibilities include new methods to study CP violation, complementary information on AM Bs 
in case Ar# s is measured first, and alternative constraints on |Vtd/Vt s |, especially for small 
values of this ratio. Finally, the theory of inclusive B decays itself can be expected to profit 
from a confrontation of the heavy quark expansion for AT Bs with experiment. In this respect 
AT Bs provides an important special case that directly probes 0{l/m\) contributions. 

As we have seen, the topic of AT Bs touches upon a rich variety of interesting physics 
issues and certainly merits the continued efforts needed to address the problems that are still 
unresolved. 
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Appendix A: Penguin Contributions to ATb s 

In the following we discuss the impact of penguin operators on the width difference Ar Ba ■ We 
will work to leading logarithmic accuracy in QCD and include the charm quark mass effects. 
For the purpose of this section we shall neglect 1 / mf, corrections, CKM suppressed modes and 
light quark masses. 

Taking gluonic penguin operators into account, the effective hamiltonian in fllUp is gener- 
alized to 
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K// = %o^£aa (64) 

V2 r= i 

where 

Ql = (biSi)v-A(CjCj) V -A Q2 = (biSj)v-A(CjCi) V - A (65) 

Q3 = (biSi)v-A(qjqj)v-A Qa = (biSj)v-A(q~jqi)v-A (66) 

Q5 = (bi.Si)v-A(qjqj)v+A Qa = (biSj) V - A (qjqi)v+A- (67) 

A summation over q = u, d, s, c is implied. Ci, . . . , C 6 are the corresponding Wilson coefficient 

functions. have already been given in fliTf). For a recent review of this subject see p2 



where further details may be found. Using our standard parameter set with // = nif, the 
numerical values are 

(Ci, . . . , C 6 ) = (-0.272, 1.120, 0.012, -0.028, 0.008, -0.035). (68) 

The calculation of the transition operator (||) using the extended operator basis is straight- 
forward and leads to 



(— = 16tt 2 B(B s -> Xev 

V r y n 



r 2 



™fo £(2) ^QCZ) 



v 7 ! - 4z • 



(2(1 - + K[ + K'{) + (1 - Az){K 2 + K' 2 + K'i) + 62^3 + K'D) [l + — )B+ 



+ (1 + 2z)(AT 2 + K' 2 + K'i -K x -K[- K'l) 



Ml 



+3 



{m b + m s 
Ml 



2 B s 

N ' 



N 



Be 



K12 are defined in ([15]) and the remaining coefficients read 



K[ = 2(iV c C 1 C 3 + C1C4 + C 2 C 3 ) K' 2 = 2C 2 C 4 
K' 3 = 2{N C C 1 C 5 + C X C 6 + C 2 C 5 + C 2 C 6 ), 



N 



(69) 



(70) 
(71) 



^ = N c Cl + iV c C 5 2 + 2C 3 C A + 2C 5 C 6 ^ = C\ + C 6 2 (72) 

^3 = 2(N C C 3 C 5 + C 3 C 6 + C4C5 + C 4 C 6 ). (73) 

These expressions represent the interference of penguin operators with the leading operators 
Qi,2 (coefficients K-) and penguin-penguin contributions (coefficients K"). Numerically they 
reduce (AT/T) Bs by 0.0114, which is about 5% of the result without penguins (AT/T) Bs = 
0.221, neglecting l/m& corrections. Note that since C3, . . . ,C$ are small, the effect of penguins 
is dominated by the while the K" are negligible. 
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Appendix B: Cabibbo-suppressed Contributions to AFb 



In this appendix we briefly consider the CKM suppressed contributions to Ar# s . They arise 
from uc (uc) or uu intermediate states in the diagram of Fig. [|. For our estimate we include 
again QCD corrections in the leading logarithmic approximation and keep charm quark mass 
effects. We neglect l/m b corrections and the small impact of penguin operators in the uu 
channel. 

The contribution from uc and uc intermediate states is then found to be 



V r J Ba ,uc g(z)r] QCD ml A c 

((2 + z)K x + (1 - z)K 2 ) (l + i-) B + (1 + 2z) (K 2 - K,) (2 - 1- ) B s 



(74) 



where Aj = V^Vi S . Compared to the leading, CKM allowed contribution with two charm 
quarks in the intermediate state, expression QT4] ) is suppressed by a factor 

2Re^ = 2\ 2 g< ±3% (75) 
A c 

Here we have used the Wolfenstein parametrization and the result that g is restricted by 
\g\ < 0.3 in the standard model [3~2|. Since the difference between ( |T4] ) and the CKM allowed 



contribution (see (|30|)) due to the different charm quark mass dependences turns out to be 
negligible numerically, relation ([T3|) determines essentially the relative importance of (|T4]) for 



(AT/T)b s - Note that the sign of ([74]) is not yet fixed because both positive and negative 
values are still allowed for g. Since g could be close to zero, the CKM suppressed contribution 
(fr"4]) might also be well below the 3% given above. In any case, it can be safely neglected. 

The contribution with two internal up quarks can be obtained from (|74]) by replacing 
2Re(A u /A c ) — > Re(A u /A c ) 2 and setting z — ► everywhere except in the argument of g(z). 
Since |Re(A u /A c ) 2 | can be estimated to be smaller than 4 • 10 -4 , the resulting expression is 
still much more suppressed than ( [74] ) and therefore completely irrelevant. 



Appendix C: Comment on r B +/rB d 

Some of the issues in the calculation of lifetime differences among B s and B^ mesons that we 
have discussed in this paper are also relevant for the prediction of tb+/tb a - We will therefore 
take the opportunity to also have a brief look at the question of the B + — Bj lifetime difference. 
In the literature this quantity has been estimated to be 0] 



, Bi - 1 + -° 5 '(200W (76) 

predicting the B + lifetime to exceed T# d by several percent. In the following we would like 
to re-examine this estimate, emphasizing the theoretical uncertainties that are involved in its 
derivation. Assuming isospin symmetry, the mechanisms that produce a difference in tb+ and 
t s d first enter at the level of dimension six operators, or equivalently at 0(l/mf), in the heavy 
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quark expansion ||. These effects are weak annihilation for the Bd and Pauli interference in 
the case of B + . As we have seen in section f|, the weak annihilation contribution to r Bd is 
very small and we shall neglect it. In this approximation the difference between t b + and r Bd 
arises only through Pauli interference and one may write 



where 



1 + 2Aix 2 B(B 



mi 
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,2 • 



The Wilson coefficients C± have been given in (|TT 

>(«) 



define the bag parameters -B} 

With m b = 4.8 GeV, m c = 1.4 GeV, A LO = 0.2 GeV and taking £[" 2 J = 1, f B = 0.2 GeV, 
one finds t b + / T Bd = 1.02, indicating a slightly longer lifetime for B + than for Bd- This 
number can however not be viewed as a very accurate prediction. In fact, the two contribu- 
tions proportional to B^ and B^ in (|77|) enter with different sign. This leads to a partial 
cancellation that has the tendency to make the result unstable. For instance, allowing the 
unphysical scale /i = 0(m b ) in the coefficients C± to vary from m b /2 to 2m b gives a range 
of 1.00 — 1.06 for the B + to Bd lifetime ratio. Switching off short distance QCD corrections 
completely (C± —* 1), the hierarchy of lifetimes would even be reversed to T B +/r Bd = 0.95, 
which is another aspect of the large sensitivity to QCD effects. An alternative way of esti- 
mating the present uncertainty is to allow a variation in the bag parameters (keeping fi = m b 
fixed). A range of B^\ = 1.0 ±0.3 is certainly conceivable, considering the uncertainties in the 
nonperturbative dynamics and from the scale and scheme dependence in the long-distance to 
short- distance matching. Assuming this, we obtain for f B = 0.2 GeV, r B +/r Bd = 1.02 ± 0.04. 
A combination of both variations, of scale and bag parameters, would even allow us to obtain 
a lifetime difference of up to 20%, r B +/r Bd ~ 1.2. Although we consider this case highly un- 
likely the point to note is that a lifetime that large could be tolerated by QCD as well as equal 
lifetimes, or even a marginally shorter lifetime for the B + . A decisive improvement of this 
situation could only be achieved by a reliable lattice calculation of B^l m conjunction with 
a next-to-leading order computation of short- distance QCD corrections to ensure a proper 
matching in renormalization scheme and scale between Wilson coefficients and hadronic ma- 
trix elements. Alternatively one could use the present measurement T B +/r Bd = 1.06 ± 0.04 
26| to constrain the bag parameters. At present such constraints appear to be of limited use, 



because of the large renormalization scale dependence of Pauli interference at leading order. 
Similar conclusions have been reached in the recent paper by Neubert and Sachrajda [j3(| . 

The authors of have modeled the bag parameters in their estimate of r B + / r Bd by 
factorizing at a low scale \ih < m, and explicitly including the leading logarithms of HQET. 
This yields 



B[ u \m b ) 



a s {m b ) 

OLs{Hh) 



-3/50 
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a s {m b ) 



-i 12/25 
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a s [m b ) 



-i 12/25 



(79) 
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Taking p h = 1 GeV this gives B[ u) (m b ) = 1.01, Bp(m b ) = 0.72 and r B +/r Bd = 1.04 (for 
fs = 0.2 GeV), favoring t b + > r Bd . However, as discussed at the end of section 0, the 
quantitative reliability of an estimate based on hybrid logarithms is not entirely clear. 
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